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RESUME. Le concept de I'invariance de jauge dans I'electrodynamique 
classique suppose tacitement que les equations de Maxwell possedent 
des solutions uniques. Mais en calculant le champ electromagnetique 
d'une particule en mouvement faisant usage de la jauge de Lorenz ainsi 
que de la jauge de Coulomb, et resolvant directement les equations des 
champs nous obtenons des solutions contradictoires. Nous concluons 
done que I'hypothese tacite de I'unicite de la solution n'est pas justifiee. 
La raison pour cette difficulte pent etre attribuee aux equations d'onde 
inhomogenes qui connectent simultanement les champs propageants et 
leurs sources. 

ABSTRACT. The concept of gauge invariance in classical electrody- 
namics assumes tacitly that Maxwell's equations have unique solutions. 
By calculating the electromagnetic field of a moving particle both in 
Lorenz and in Coulomb gauge and directly from the field equations we 
obtain, however, contradicting solutions. We conclude that the tacit 
assumption of uniqueness is not justified. The reason for this failure 
is traced back to the inhomogeneous wave equations which connect the 
propagating fields and their sources at the same time. 
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1 Introduction 

"The principle of gauge invariance plays a key role in the standard model, 
which describes electroweak and strong interactions of elementary parti- 
cles." This statement is quoted from an article by J. D. Jackson and L. B. 
Okun [1] which emphasizes the importance of the principle and delineates 
its historical evolution. The concept emerged from classical electrody- 
namics where the electromagnetic field expressed in terms of potentials: 
E = —V(p — dA/cdt, B = rot A does not change under the transforma- 
tion: A + Vx, (j) ^ (j> — dx/c dt. Since div A — > div A + Ax and x 
is an arbitrary function, the divergence of the vector potential can seem- 
ingly be chosen arbitrarily without influencing the fields. This feature 
was exploited to de-couple Maxwell's inhomogeneous equations by either 
choosing div A + dip/cdt = (Lorenz'' gauge) or div A = (Coulomb 
gauge) . The solution for the fields should be entirely independent of this 
choice. 

There is. however, a tacit assumption behind the formalism of elec- 
trodynamic gauge invariance: Maxwell's equations must have unique so- 
lutions, otherwise it is meaningless to talk about potentials from which 
fields may be derived. In reference [1] it is said: "It took almost a cen- 
tury to formulate this nonuniqueness of potentials that exists despite 
the uniqueness of the electromagnetic fields." To our knowledge it was 
never attempted to prove that the electromagnetic field resulting from 
a solution of Maxwell's equations is actually unique, it was just taken 
for granted. If there were no unique solution of Maxwell's linear system 
of first order field equations, gauge transformations on (nonexisting) po- 
tentials would be irrelevant. 

In this paper we start with the usual assumption that unique solu- 
tions of Maxwell's equations do exist and try to calculate them both 
with the help of potentials (Sections 2-4) and directly from the field 
equations (Section 5). For the electromagnetic field of a moving parti- 
cle we find, however, contradic-tiiig solutions depending on the method 
used. In Section 6 we show that the standard Lienard-Wiechert fields, 
which satisfy Maxwell's source-free equations, cannot be considered as a 
unique solution of Maxwell's inhomogeneous equations. Thus, we infer 
that the tacit assumption concerning the existence of unique solutions is 
not justified in general. 

^In [1] it is pointed out that Ludwig Valentin Lorenz published more than 25 years 
before Hendrik Antoon Lorentz what is commonly known as the "Lorentz condition" . 
To give proper credit to Lorenz we use in this paper the term "Lorenz gauge" . 
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The reason for this failure is discussed in Section 8 where we come to 
the conclusion that the mixture of elliptic and hyperbolic equations, as 
formulated by Maxwell, does not permit a physical solution for moving 
point sources. 



2 The electromagnetic field of a moving particle calculated in 
Lorenz gauge 

The electromagnetic field produced by a moving particle is calculated [2] 
from Maxwell's equations: 

divE = 47rp (1) 



rotE = --- (2) 



divB = (3) 



47r ^ 1 dE 
rotB=-pv+-- (4) 

where p is the charge density located in a narrow region round a center 
moving with velocity v. The total charge of the particle is the integral 
over all space: 

JJJp {x', t) d^x' = e (5) 

The usual method to solve equations (1-4) is to adopt the potential 
ansatz: ^ 

1 3 A 

E = — - V(f) , B = rot A (6) 

c ot 

which satisfies automatically equations (2) and (3) and leads to two 
second order differential equations: 

1 dA 

A(t)+ -div— = -4 TT p (7) 
c ot 

1 d"^! 47r ^ _ / . 1 
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The freedom in the choice of the divergence of the vector potential 
is exploited to de-couple equations (7) and (8). Adopting the Lorenz 
condition: 

divA+-^ = (9) 
c ot 

one obtains two wave equations of the same structure: 

^ 1 47r ^ 

The solution of (10) subject to the boundary condition that the scalar 
potential vanishes at infinity is: 

</.(x,t) = yy y ^^^''^~^^ dv , r=\x-x'\ (12) 

where the position of the charge density is to be taken at the retarded 
time: 

t' = t-- (13) 

c 

The advanced solution with t' = t + R/cis excluded on physical grounds. 

The integral (12) may be carried out by employing the 5 - function 
formalism: 

+00 

— 00 

Substituting this into (12) yields: 

+00 

cP{x,t)= J j j j pS^^5 [t' -t + ^ Sx'dt' (15) 

— 00 

The integration over all space for a point-like particle results with (5) 
in: 

+00 

{x,t) = j ^5 -t + ^ dt' (16) 
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where R expresses now the distance between the field point and the 
position of the charge x' {t') at the retarded time. Changing to the 
variable 



u = t'-t+ 



R {x, x' {t')) 



du ^ 

Ij^ ~ cR dt' 



R dx' it') dx'it') 



dt' 



R 



we may integrate (16) and obtain the result: 

(j) {x,t) = 

Similarly, we find from (11): 
A{x,t) = 



v{t') 
(17) 



(18) 



iR 



(l-fl) 



(19) 



Solutions (18) and (19) are the well-known retarded Lienard-Wiechert 
potentials. With the differentiation rules resulting from (13): 



d 



1 



d 



1 d 



d 



dt 1 _R:v dt' A dt' ' dx 

Rc 

one obtains with (6) for the fields: 



5 (/ (i')) 



df 1 dR dR _ R 
dt' c dx ' dx XR 
(20) 



E {x, t) = e 



B{x,t)- 



1_ / _ _ 



R^ 



dv 



dt' ) c^R\^ dt' 



R X V 



1 - dv\ 



1-^ +-R.— \ + 



R 



dv 



dt' J c^R'^X'^ dt' 



(21) 



where t' = t — ^. 



3 Solution in Coulomb gauge 

The fields as given by (21) must be the same when they are calculated 
in Coulomb gauge by adopting the condition: 

divA = (22) 
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Equations (7) and (8) become: 

= - 4 TT p 



(23) 



, J. I d'^A 47r ^ 1 o> _ , 

Solution of (23) yields the instantaneous Coulomb potential: 

(j>c {x,t) = - , r=\x-x'{t)\ (25) 
r 

which substituted into (24) results in a wave equation for the 'Coulomb' 
vector potential: 



1 d'^Ac 4 7r 



d 



(26) 



The vector r = x — x' (t) denotes now the simultaneous distance between 
charge and field point. The first term on the right-hand-side of (26) yields 
a contribution to the vector potential which is identical with (19): 



Aci {x,t) 



(27) 



The solution of the wave equation for the second part of the vector 
potential may be written in the form: 



+ 00 



Ac2 {x,t) 



47rc 
f=s-x'{t') 



TTTT ^ 5{t' -t+- \ —dt' 



dt 



R 



R 



(28) 



when we employ & 5 - function as in the previous Section. Here wc have 
used the integration variable s in distinction of the position x ' (t) of the 
charge. The contribution of Ac2 to the electric field is with (6): 



dt'^ 
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and the contribution to the magnetic field may be written as; 



B2 {x,t)- 



Ait c 



V (^) (4) Ssdt' 



R 



dt 



Attc 



5' 



cE? E? 



dt' 



d^sdt'{30) 



Performing a partial integration over t' we obtain the solution in the 
form: 



B2 {x,t) = 



47rc 



R 
RP 



d_ 



R 
c 



(31) 

Since (27) yields already the magnetic field as given by (21), the 
contribution (31) must vanish which is not likely to occur: The finite 
cross-product is to be integrated with different weight, so that both 
terms in (31) cannot vanish simultaneously. Furthermore, the condition 
for identical electric fields in Coulomb and in Lorenz gauge: 



1 dAc2 
c dt 



(32) 



where (pLW denotes the potential (18), cannot be satisfied, if (28) is not 
an irrotational field. Even if the second term in (32) could be written 
as a gradient: dAc2/cdt = V<?!)2, the condition: (pLW = 4>C + 4'2 + const 
would still be violated, since dAc2/dt, and thereby ^2, would depend 
on the acceleration according to (28), which is not the case for ^lw nor 
<l>c- 

In order to quantify these qualitative considerations the integral (28): 



f=s-x'{t') , 



3 r {r • v) 



V 



R = s — X 



d^ 

_ R R 



(33) 



may be evaluated analytically for the case of a constant velocity of the 
charge: 

x' {t') = XQ + vot' = xo+ vo (^t - (34) 
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The integration variable s may be replaced by R so that the vector r 
may be written as: 



-» R 

f = R + X - (xq + Vot) -\ Vo 

c 



(35) 



We assume that the charge moves along the z-axis of a coordinate system 
having its origin at _R = 0. The z-component of the vector potential 
evaluated on the z-axis becomes then; 



Ac2z = 



4tt c 



1 3 {R,+PR + zY 



f_R 

R ' ' c 



r^^R^ (1 + /3^) +2(3RR^ + 2z {PR + R^) + 



(36) 



where z {t) denotes the distance between the field point and the position 
of the charge at time t. Using spherical coordinates: 



R = R sm.9 cos i+R sin 9 sin ip j+R cos 9 k , S R = R? sm.9 dR d9 dip 

(37) 

expression (36) becomes: 

2 TT TT CO 


„2 d2 /i , o2 



1 3 {R {f3 + cos9) + zy 



R sm9dRd9dip 



R^ (l + /3^ + 2/3 cos6') + 2Rz {(3 + cost 



(38) 



The integration over ip yields a factor of 2 vr, since the intc^grand is inde- 
pendent of ip. Upon indefinite integration over R and 9 one obtains: 



Ac2z 



(39) 



i?2 {R+{(3R + z) cos9) 



'2'CZ {f3R + zy ^JR^ (l-F/?2 + 2/3cos6')-h2i?z (/? cos 6*) + 2^ 

The integral vanishes both at i? = and at i? = cx). It is singular 
at i? = 2;/ (1 — /3) , 6 = TT. Close to the singularity we expand it by 
substituting R = z/ {1 — (3) — ze\ , 6 = n — €2 and obtain in lowest 
order: 



2cz Ve2 4 62 



(40) 
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Obviously, the integral assumes no definite value when we go to the 
limits ei = £2 = , as it does not converge absolutely. 

If wc perform the same calculation on (29) we obtain in addition to 
undefined terms a diverging contribution: 

evl l-f3 
2 ei 

which is also encountered when we calculate (31). 

Prom these results we must conclude that Maxwell's equations do not 

yield a physical solution for the fields of a moving particle in Coulomb 
gauge. Furthermore, the undefined fields as derived from (25), (27) and 
(28) by using (6) disagree definitely with the fields as given by (21) in 
Lorcnz gauge which arc well defined. Similar conclusions were reached 
by Onoochin [3] without evaluating the integral (33) explicitly. 

Jackson has attempted [4] to find an approximate 'quasistatic' so- 
lution in Coulomb gauge which should be valid for velocities v <^ c m 
a region very close to the particle where retardation may be neglected. 
We discuss this attempt in Appendix A and show that it also leads to 
an inconsistency. 

In a recent paper [6] Hnizdo has given a solution of (26) which is 
based on the gauge function as defined in equation (66) below in Section 
7, and on the formal solution (67) for the gauge function. The second 
formal solution (70) is ignored in this consideration. Hnizdo arrives at 
a similar result as we found in (40), but he establishes uniqueness by 
applying a regularization procedure which can hardly be justified from 
a mathematical point of view. We discuss his approach in Appendix B. 



4 Helmholtz's ansatz 

Having obtained contradicting solutions in Lorenz gauge and in Coulomb 

gauge we infer that Maxwell's equations contain an inconsistency which 
does not permit to find a unique solution for the fields. In order to 
trace this problem, we employ Helmholtz's theorem which states that 
any vector field may be expressed as the s\im of a rotational and an 
irrotational field. This was already used in the ansatz (6). Now we 
apply it to the electric field: 

E = rotU -W(f> , divU = (42) 
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Substituting this into (2) we obtain: 



AU = 



1 dB 
c 'dt 



(43) 



Taking the rotation of (4) and inserting (42) yields a second Poisson 
equation: 



Its solution for a point charge is: 

1 dU_ 
c 'dt 



47r „ 

Vp XV 

c 



B 



(44) 



(45) 



Substituting this into (43) yields a wave equation for the vector potential 
of the electric field: 



e_ d_ 

? dt 



V X r 



(46) 



which has the retarded solution: 
U 



d_ 

df 



47rc2 

r = s — x'{t') , R — s — x 
For the magnetic field wc obtain with (45): 



'=t--S 



V X r 



d-^s e 



V X r 



(47) 



(48) 



and the electric field as derived from (47) with (42) becomes: 



iTTC^/// 



R 

c 



V X r 



+ 



d_ 
dt' 



V X r 



t'=t- 



e r 



(49) 



where we have added the Coulomb field which results from insertion of 
(42) into (1). We note that neither of the expressions (48) and (49) agrees 
with the fields as calculated in Sections 2 and 3, because the fields derived 
from the Hclmholtz ansatz (42) depend on the second time derivative of 
the velocity. Assuming a constant velocity of the particle one could also 
show that the integrals (48) and (49) actually diverge. 
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5 Direct solution of the field equations 

The two types of potential ansatz (6) and (42) resulted in different so- 
lutions for the fields. We, therefore, want to calculate the fields directly 
from (1-4) without using any potential ansatz. By elimination of the 
electric and the magnetic field, respectively, we find the wave equations; 



AS 



c 



1 d^B 
"2 -g^ 



An 



(50) 



-5 1 d^E An d , „ 



(51) 



df^ C2 dt 

By applying the standard method of solving this kind of wave equation, 
as it was described in Section 2, it can be shown that the ensuing solution 
of (50) and (51) yields exactly the fields as given by (21). However, by 
deriving the hyperbolic equations (50) and (51) we have ignored the fact 
that Maxwell's equations are actually a mixture of hyperbolic and elliptic 
equations which became very obvious in the previous Section. In order 
to take this into account we split the electric field into its irrotational 
and its rotational part: 

E = Eg + Er (52) 

The rotational part does not enter equation (1). The irrotational part is 
just the quasistatic Coulomb field which does not propagate: 



E„ 



e r 



, r = x-x'{t) 



(53) 



The rotational part obeys the wave equation: 



/^Er 



which has the retarded solution: 



dt 



(pv) 



C2 5*2 



(54) 



Er {x,t) = 



1 d {p is, t') V it')) 

C2 dt' 



47rc2 5f'2 



(55) 



|£-£| IX — S\ 



Adding (53) one obtains the electric field as it was derived in Section 
3 in Coulomb gauge, whereas the retarded solution of (50) yields the 
magnetic field as it was derived in Section 2 in Lorenz gauge. It was 
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shown in Section 3 that the integral (55) diverges. Hence, it does not 
represent a physical solution for the rotational part of the electric field. 

This analysis shows that the inconsistency inherent to Maxwell's 
equations is not an artefact produced by employing a potential ansatz. 
It seems to result from the mixture of hyperbolic and elliptic differential 
equations for the fields, as they were formulated by Maxwell. Only in 
Lorenz gauge the elliptic equations are completely removed so that there 
is seemingly agreement between the solutions of the hyperbolic equations 
(50) and (51) and the hyperbolic potential equations (10) and (11). 



6 The inhomogeneous wave equations for a moving point 
source 

The discrepancies encountered in Sections 3-5 are apparently related 

to the fact that Maxwell's set of equations mixes hyperbolic and elliptic 
structures so that unique solutions are not possible. In view of this find- 
ing it is somewhat surprising that in Lorenz gauge the elliptic features 
seem to be removed altogether so that equations (10) and (11) yield the 
unique solutions (18) and (19), provided the advanced solutions are sup- 
pressed on physical grounds. A closer look at the inhomogeneous wave 
equation (10) reveals, however, that the elliptic character is still there, 
but concealed in the inhomogeneity. If it is brought out, it turns out 
that the solution (18) cannot be considered as unique. 

In order to sec this we employ a different method of solution than 
that used in Section 2. Due to the linearity of (10) one may split the 
potential into two contributions: 

(j> = (j>o + (j)^ (56) 

The wave equation may then be split into a Poisson equation: 

A0o=-47rp (57) 

and into another wave equation: 

1 a2</.i 1 d^,i>o 



where the solution of the elliptic equation (57) enters as an extended 
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source'' Adding (57) and (58) the wave equation (10) is recovered. The 
retarded solution of (58) is: 



01 {X, t) = 



t'=t--^ 



R 



r=s ■ 



x'{t') 



R 



(59) 



where we have substituted the instantaneous Coulomb potential (25) 
resulting from a solution of (57). Carrying out the differentiation in (59) 
we have: 



01 {x, t) 



47rc2 



1 



— r 



dv 



3 (r-v) 



(60) 

This integral depends on the acceleration which is not the case for the 
potential (18). At constant velocity the integral is very similar to (33) 
and we know from the calculation in Section 3 that it has no defined 
limiting value according to (40). 

Similar considerations apply to the inhomogcncous wave equation 
(11), the solution of which may formally bo written as: 



A (x, t) = — 

cr 



d^s 



by applying the same method of splitting the vector potential into two 
parts. If result (61) is substituted into (6), the fields would depend on 
the third derivative of the velocity, which is not the case according to 
(21), so that (61) is incompatible with (19). 

We must conclude then that the potentials (56) and (61) disagree 
with the Licnarci-Wiechert potentials (18) and (19) which turn out not 
to be a unique solution of the inhomogeneous wave equations (10) and 
(11), even if the advanced solutions are suppressed. As a matter of fact, 
equations (10) and (11) have no physical solution judged from our results 
(60) and (61) which do not admit a defined limiting value. 

There is a direct way of showing that the Lienard-Wiechert result, 
which leads to the fields (21), satisfies only Maxwell's homogeneous equa- 
tions, but the inhomogeneities are not taken into account properly. Let 

^The ansatz (56) was also used in [2], but it was erroneously assumed that 
d'^(t>o/dt'^ in (58) may be neglected. 
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us consider Green's first identity: 

JJ<l>V(f)-(fx = j j j ((/)A(?i+ |V^|^) d^x (62) 

The surface integral on the left-hand-side vanishes at infinity both for 
the Coulomb and for the Lorenz potential so that the volume integral 
over all space on the right-hand-side must vanish as well. Substituting 
the Coulomb potential (25) together with (23) one obtains an integral 
equation which must be satisfied by the charge density: 



oo 

Aire j p (r) -|- rfr = 



(63) 



where spherical coordinates centered around the position of the charge 
were used. If one inserts the Lienard-Wiechert potential (18) together 
with (10) into (62), one obtains: 

oo TT 



(64) 

This integral equation depends now on the velocity which may be easily 
verified by choosing a constant velocity so that (18) yields: 

<I>LW = r , r = x-{xo + vt) (65) 

(v ■ r/cf + (1 — v'^/c^) 

Both integral equations (63) and (64) cannot be satisfied by the same 
function p (r), unless the 'shape' of the point charge would depend on 
the velocity, as suggested by Onoochin in Reference [3]. 



7 Treinsformation of the Lorenz potentials into Coulomb po- 
tentials 

From the results obtained in Sections 2 - 6 it should be evident by now 
that a unique gauge transformation, which transforms the Lorenz poten- 
tials of a point source into the corresponding Coulomb potentials, cannot 
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exist. We finally want to show this explicitly. The gauge transformation 
is effected by a generating function 



<i)c = 4>L- 



1 dx {x, t) 



Ac = Al + Vx {x, t) 



c dt 

By integrating the first relation over time one obtains immediately: 

t 

X{x,t) = j c {(f>L - (t>c) dt + xo (x) 

to 



(66) 



(67) 



where the Lorenz potential is given by (18) and the Coulomb potential 
by (25). The gauge function must also satisfy the Poisson equation: 



Ax {x, t) = div Ac — div Al 



(68) 



which follows from the second relation in (66). Equation (68) has the 
instantaneous solution: 



X (x, t) 



" d(j>L {s, t) d^s 



dt 



\s — x\ 



(69) 



where we have substituted the conditions (9) and (22) into (68) . If we 
insert the Lorenz potential as given by (18) and apply the first differen- 
tiation rule in (20), we find: 



X {x, t) 



An c 



A3 i?3 



c 



R- 



dv 
dt' 



d^i 



R = g-S'{t') , t' = t- 



■x\ 



\s-x\ 
(70) 



As this expression depends not only on the velocity, but also on the 
acceleration, it is not compatible with expression (67). Furthermore, the 
integral (70) has no unique limiting value, but depends on the chosen 
integration variables. In order to see this we assume a constant velocity. 
The Lienard-Wiechert potential (18) becomes in this case: 



(l>L {X, t) = 



{v ■ f/c)'^ + (1 ~ v'^/c^) 



f^x-x'it) (71) 
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Substitution into (69) yields: 



V ■ s 



(72) 



{v ■ s/cf + (1 — v'^/c^) s' 



,2li \r-s\ 



where we have chosen a coordinate system with its origin at the position 

x' {t) of the charge. If we change to the integration variable s' = s — r, 
which is equivalent to shifting the origin of the coordinate system to the 
field point x, we obtain instead: 



Apart from a common logarithmic singularity at infinity the integrals 
(72) and (73) are conditionally convergent and assume diff'erent limiting 
values which was verified by calculating them in spherical coordinates. 
Numerical calculations of (72) and (73) in cylindrical coordinates yield 
still diflferent limiting values. None of these results agrees with (67). 

In accordance with the previous conclusions reached in this paper we 
infer that no unique function x exists which would transform the Lorenz 
potentials of a point source into the corresponding Coulomb potentials. 
Hence, the principle of 'gauge invariance' is not applicable to classical 
electrodynamics. 



8 Discussion 

The nature of the inconsistencies encountered in Sections 2 - 5 is appar- 
ently connected to the feature of Maxwell's equations of mixing elliptic 
and hyperbolic structures. Even if the equations are reduced to wave 
equations like (10) and (11) or (50) and (51), the elliptic character is 
still there in form of the inhomogeneity and may be made visible by the 
method of solution employed in Section 6. There we were compelled 
to conclude that the inhomogeneous wave equation does not have, as a 
matter of fact, a unique solution, or even leads to unphysical diverging 
solutions, at least in case of a moving point source. In principle, it is well 
known that the inhomogeneous wave equation has an infinite manifold 



X {x, t) = 



(73) 




V ■ {s' + r) 



{v- (s ' + r) jcf + (1 - t>Vc2) (s'2 j^Tfi-2s> -r) 
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of solutions, but it is generally believed that suppression of the advanced 
solutions reduces it to a physical solution, the properties of which are 
uniquely determined by the behaviour of the source. According to our 
analysis in Section 6 we must maintain, however, that the inhomoge- 
neous wave equations do not correctly describe the physical process of 
wave production by a moving point source. 

Our result is not too surprising, if we realize that the inhomoge- 
neous wave equations (50) and (51) relate the measurable fields at a 
certain location with temporal changes in a remote source at the same 
time. Although both the fields and the sources in equations (1 - 4) were 
diff'erentiated at time t when, e.g., the electric field was eliminated to 
obtain (50), the retarded solutions (21) require that the differentiation 
of the sources is dated back to the earlier time t — R/c. However, the 
source may be an extinguished star the light of which we see only now 
at time t. It makes no sense to differentiate a source not existing any 
more which has no influence whatsoever on the light we see after a bil- 
lion years. When we differentiate Maxwell's equations now to obtain 
the wave equations for the fields, we treat the temporal changes in the 
sources as if they would happen now at time t. In the retarded solu- 
tions, however, we date back the change in the sources to a remote past. 
This procedure is inconsistent, but inescapable due to the structure of 
Maxwell's equations. 

If the same procedure would be applied to acoustic waves, one would 
encounter similar inconsistencies. Instead, from the hydrodynamic equa- 
tions one derives linearized homogeneous wave equations for the pressure 
and the fluid velocity. These are solved by imposing suitable boundary 
conditions which are determined, e.g., by the oscillating membrane of a 
loudspeaker. Maxwell's equations, however, lead to inhomogeneous wave 
equations which connect the travelling fields with the source at the same 
time, a contradictio in adjecto. This became quite obvious in Section 
5 where equation (51) predicts that the total electric field is produced 
in a point-like region and travels within a finite time to the field point 
where an observer may be placed. On the other hand, equation (53) 
predicts that part of the field has already arrived there instantaneously, 
as soon as any change in the source occurred. This inconsistency cannot 
be resolved without altering equations (1 - 4). 

It is well possible that Maxwell was fully aware of this problem, 
because in his 'Treatise' [7] he did not derive an inhomogeneous wave 
equation. He used the Coulomb gauge and derived equation (24). Then 
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he argued that in the 'ether' there does not exist any current or free 
charge. This way he was left with: 

, J. 1 1 , , „ 

Now he committed a formal inconsistency by concluding that the van- 
ishing of the Laplacian of the scalar potential justifies to omit also the 
gradient of the potential. This is, of course, not true close to a charged 
body. In other words, he omitted part of the 'displacement current', 
which was invented by him in the first place, in order to 'derive' a ho- 
mogeneous wave equation for the vector potential. At last he suggested 
that this equation should be solved by imposing Cauchy-type boundary 
conditions on A and dA/dt. The result would be a travelling vector 
wave from which the measurable fields could be derived with (6) at 
any place and time where the wave has arrived. In his last Chapter 
XXIII (Article 862) Maxwell discusses Riemann's inhomogcncous wave 
equation which is formally the same as Lorenz's equation (10). Taking 
reference to Clausius he states that Riemann's formula is not in agree- 
ment with "the known laws of electrodynamics" . Maxwell's method of 
using only a homogeneous wave equation to describe electromagnetic 
waves is still in practical use, when the radiation emitted by an antenna 
is calculated. Only the homogeneous wave equation is used together 
with plausible boundary conditions resulting from a physics which goes 
beyond Maxwell's equations. 

In parentheses we remark that Einstein gave his famous paper of 
1905 the title: "Zur Elcktrodynamik bewegter Korper" , but the 'moving 
bodies' which carry charges and currents are not treated in his analysis. 
He deals only with Maxwell's homogeneous equations which do not lead 
to contradictions. The instantaneous Coulomb potential is left out of 
the consideration. 

For slowly varying fields Maxwell's equations do make sense. They 
describe correctly the phenomenon of induction in transformers where 
only the instantaneous fields come into play and where the displacement 
current is negligible. When a condenser is charged up, the displace- 
ment current must be allowed for, but then it is only the instantaneous 
Coulomb field which matters in practice. A quasistatic 'near field' the- 
ory can be carried through satisfactorily, but amalgamating it with wave 
phenomena, in the way Maxwell has tried it, leads to the contradictions 
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which we have demonstrated. These were also recognized by Dunn [8], 
but not worked out in detail. 



Appendix A 

The wave equation (26) for the second part of the 'Coulomb' vector 
potential: 



may be taken as a Poisson equation which has the formal solution: 

Ac2 = 



Attc 



V 3 f {v ■ r)^ (Ps 



Close to the charge and at small velocity t; <C c the second term may 
be expected to be negligibly small so that the first integral could be 
considered as an approximate solution of (Al) with limited applicability. 

This attempt to obtain a 'quasistatic' solution was pursued by Jack- 
son [4] in order to find the interaction Lagrangian between two particles 
moving at nonrelativistic velocities. He chose a coordinate system cen- 
tered at ^' = and performed a partial integration: 

V ■ s\ d f 1 d f f f V ■ s d^s 



c3 



ds \ \s — x\ I dx 



Now the integration was straightforward and yielded: 



r e d f V ■ x\ el v x (v ■ x) \ , , 



In order to obtain the total vector potential, the unretarded contribution 
from expression (27) must be added: 

T- e I V X (v- x)\ / . ,.x 

^c = ^{t^ + ^-ir^ (A.5) 



2 c V If I 



X 
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It turns out, however, that the solution (A. 4) does not satisfy the 
Poisson equation from which it was calculated. Substituting (A. 4) into 
the left-hand-side of (A.l) and ignoring the second time derivative yields: 



The distance vector x pointing from the origin to the field point may 

be identified with r, since the position of the charge was assumed at the 
origin. The second term of (A.6) equals the right-hand-side of (A.l), but 
the first term yields a ^ - function: 



which remains unaccounted for in (A.l). 

The reason for the discrepancy is that the first integral in (A. 2) does 

not converge absolutely, as wc have shown in Section 3. Consequently, 
the operations of partial integration as well as interchanging the sequence 
of differentiation and integration in (A. 3) are not permitted and lead to 
an incorrect result. The first integral in (A. 2) has, in fact, no defined 
limiting value as is obvious from expression (40). 

In a private communication Professor Jackson explained how he 
could obtain (A5) from an expansion procedure applied on the Lienard- 

Wiechert potential - which is based on the Lorenz gauge - and that this 
was actually being done by Darwin who derived (A. 5) in 1920 [5]. 



Appendix B 

Hnizdo's article Potentials of a uniformly moving point charge in the 
Coulomb gauge [6] was apparently motivated by Onoochin's objections 
against 'mainstream' electrodynamics as published in Ref. [3] . Onoochin 
reached conchisions which are similar to those arrived at in the present 
paper, namely that the electrodynamic field depends on the choice of the 
gauge. He speculates that the 'shape' of the electron could depend on its 
velocity, an idea which was already pursued (unsuccessfully) by Lorcntz. 
Hnizdo tries to resolve the problem by using a certain regularization 
procedure when the integral (33) is evaluated. 




(A.6) 




(A.7) 
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Both authors do not emphasize that this type of integral is condition- 
ally convergent. This property may lead to the known fact that the value 
of the integral depends on the sequence of integration, on the chosen co- 
ordinate system, or on the position of the origin. A typical example is 
given in Bronstein-Semendjajew's Taschenbuch der Mathematik, Verlag 
Harri Deutsch, Frankfurt, on page 347: 



The result depends on whether one integrates first over x and then over 
y, or vice versa. Although one obtains a finite result in both cases, its 
value is not unique. 

Hnizdo uses the formal solution (67) to obtain his formulas (16 - 
18) for the difference between the vector potentials in Coulomb and in 
Lorenz gauge. He claims that these formulae do not have singularities, 
but this is actually not true. The y-component for the difference is, e.g.; 

Acy - A,, = ^^ij-^ (Vc - V,) (B.2) 

If this expression is expanded around the point 
{x = vt + e\, y = -|- e2 , ^ = -|- 63), one obtains: 

Acy - Aj^y = {Vc - Vl) (B.3) 

This result is similar to our result (40) where we concluded that no 
limiting value exists. In fact, (B.3) can assume any value between zero 
and infinity depending on the way how one approaches the limits ei 
0, €2 ^ 0, £3 ^ 0. 

Hnizdo avoids the ambiguity by using a regularization procedure as 
defined in his equation (32). It amounts to using spherical coordinates 
centered around the charge point. This prescription is, however, arbi- 
trary and lacks a rigorous justification. If one deviates from it, one finds 
non-uniqueness as demonstrated in our comparison of equations (72) and 
(73). It should be noted that Hnizdo's regularization would not work, 
when applied to the expression (70) for the gauge function which is in- 
compatible with (67). This demonstrates again that Maxwell's equations 
do not have a unique solution (for moving point charges) representing a 
measurable well defined field. 
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